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The dependence on temperature of the spin magnetization of a two-dimensional electron gas at 
filling factor unity is studied. Using classical Monte Carlo simulations we analyze the effect that the 
equivalence between topological and electrical charge has on the the behavior of the magnetization. 
We find that at intermediate temperatures the spin polarization increases in a thirty per cent due 
to the Hartree interaction between charge fluctuations. 
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The physics of a two-dimensional electron gas (2DEG) 
in a magnetic field ,B, is determined by the number of 
Landau levels occupied by the electrons. Since the degen- 
eracy of the Landau levels increases with B, all the elec- 
trons can be accommodated in the lowest Landau level 
for strong enough fields. At filling factor unity, v = 1, 
and zero temperature, T = 0, the ground state of the 
2DEG is an itinerant ferromagnet. |QJ|] The Zeeman cou- 
pling between the electron spin and the magnetic field de- 
termines the orientation of the ferromagnet polarization. 
For zero Zeeman coupling the interaction between the 
carriers produces a spontaneous spin magnetic moment. 

Recently the magnetization of the quantum Hall sys- 
tem was measured || for filling factors 0.66 < v < 1.76 
and temperatures 1.55K < T < 20K. At v = 1 and very 
low temperatures the system is fully polarized, while for 
other filling factors the magnetization is reduced. The 
demagnetization for v ^ 1 is related to the existence of 
Skyrmions in the system. [Q Further experiments have 
verified the existence of Skyrmions using transport, ca- 
pacity and optical experiments. |5||| 

For v = 1, the temperature dependence of the mag- 
netization M(T), has been measured using NMR tech- 
niques ||] and magneto-optical absorption experiments. 
H Different theoretical approaches have been used for 
the study of M{T): i) Read and Sachdev have stud- 
ied the N — ► oo limit of a quantum continuum field the- 
ory model for the spin vector field, m(r). This model 
describes the long-wavelength collective behavior of the 
electronic spins. This work has been extended by Timm 
et al. H The field theory is expected to be accurate at low 
temperatures and weak Zeeman coupling. Using SU (N) 
and O(N) symmetries in the large N limit, and using 
the spin stiffness p, as a parameter Read and Sachdev 
obtained results for M(T) which are in reasonable good 
agreement with the experimental data at low temper- 
atures. ii)Kasner and MacDonald || calculated M(T) 
using many-particle diagrammatic techniques which in- 
clude spin- wave excitations and electron spin- wave inter- 
action. This theory is a good improvement on the one- 
particle Hartree-Fock theory, but it gives a polarization 
for the system too high compared with the experimen- 



tal one. Progress in the diagrammatic approach, includ- 
ing temperature dependence screening, has been done by 
Haussman. |Io|| iii) The dependence of the magnetization 
on T has been also obtained by exact diagonalization 
of the many-particle Hamiltonian for a small (up to 9) 
number of electrons on a sphere. jll],[l2) These calcula- 
tion have important finite size effects at low temperature 
and weak Zeeman coupling, iv) Finally, quantum Monte 
Carlo (MC) techniques have been used in order to cal- 
culate M(T) for a spin 1/2 quantum Heisenberg model 
on a square lattice with exchange interactions adjusted 
to reproduce the spin stiffness of the quantum Hall fer- 
romagnet. jl3| These calculations are essentially exacts 
and probably are free of finite size effects. 

A unique property of the quantum Hall ferromagnets is 
the equivalence between the topological charge associated 
with m and the electrical charge, g] This equivalence 
make the Skyrmions to be the relevant charged excita- 
tion of the 2DEG at v = 1. Charge conservation implies 
that at a given v the integral of the topological charge 
over all the space should be constant, independently of 
T. At v — 1 this constant is zero. A spin vector field 
texture produces a modulation of the topological charge 
density. Spatial spin fluctuations increase with tempera- 
ture and produce a modulation of the topological charge 
density. In this way thermal fluctuations can produce 
a strong charge fluctuation. Because of the equivalence 
between topological and real charge, the modulation of 
the charge density costs Hartree energy. The models de- 
scribed above: diagrammatic techniques, quantum field 
theory and quantum MC calculations, do not take into 
account the Hartree contribution to the energy of the 
ferromagnet. 

In this work we study the effect that the Hartree en- 
ergy has on the temperature dependence of the magne- 
tization. We perform classical MC simulations of M(T) 
for the energy functional of the Hall ferromagnet. From 
the comparison of the results with and without Hartree 
energy we conclude that at moderate temperatures the 
inclusion of this term modifies the value of the magne- 
tization up to a thirty per cent. For realistic values of 
the Zeeman coupling, we find that at intermediate-high 
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temperature the Hartree energy is near one third of the 
Zeeman and exchange energies. 

The classical model can not describe correctly the low 
temperature behavior of M(T). The classical dynamics 
of the electron spins neglects several effects, in particular 
the quantum description of the spin-density waves which 
is extremely important for describing M(T) at low tem- 
peratures. The inclusion, in the classical model of a tem- 
perature dependent low-energy cutoff simulates quantum 
effects. As we show latter, in the Hall ferromagnet this 
cutoff appears naturally when discretizing the continuum 
model. In this work we are interested in the effect that 
the inclusion of the Hartree term has on M(T). We ex- 
pect that in the quantum Heisenberg model, the Hartree 
energy should have the same effect than in the classical 
model. 

The long-wavelength and low-energy properties of the 
v = 1 Hall ferromagnet can be described by a functional 
E of the unit vector field m(r) which describes the local 
orientation of the spin magnetic order. The functional 
E has three terms the gradient leading or exchange 
term E Xl the Zeeman term E z and the Hartree term E c . 
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here t = g* [IbB /2 is the Zeeman coupling strength, £ is 
the magnetic length, v(|r— r'|) is the Coulomb interaction 
and n(r) is the charge density and it is given by 

ra(r) = ^-e uli m(r) • [d u m(r) x 9^m(r)] . (5) 

In this continuum model the total topological charge, Q, 
is given by the integral over all the space of n(r), and it 
represents the number of times m(r) winds around the 
sphere S 2 . Skyrmions are non trivial extrema solutions 
of the functional E with Q ^ 0. The continuum model 
gives a good description of the Skyrmions with moderate 
and large spins. 

In order to perform classical MC calculations we dis- 
cretize the space using a square lattice with a lattice pa- 
rameter a^. The lattice parameter corresponding to an 
electron per unit cell is a x = ^/2tt£. The energy func- 
tional has the form, 

a 2 1 

<i,3> 1 I i,j 

(6) 

Here, fij is the unit vector at site i, qi is the topological 
charge attached to the unit cell i and Vij is the Coulomb 



interaction between two charges unity distributed uni- 
formly in the cells i and j. The cell i is defined by the 
points l:(i x ,i y ), 2:(i x +1 ,i y ), 3:{i x +1,^+1) and 4:{i x ,i y 
+1). The expression of as a function of the unit vectors 
at the points 1-4 is |l6| 

gi = ^-{(tTA)(ni > n ai n 3 ) + (crA)(n 1 ,na 1 n4)} , (7) 

where (aA)(Sli, ^3) denoted the signed area of the 
spherical triangle with corners Oi , O2 and Q,^. Apart 
from a set of 'exceptional' configurations of measure zero, 
this prescription for the topological charge yields well de- 
fined integer values for the total topological charge. For 
smooth spin texture the continuum and discrete expres- 
sion for the density of topological charge give the same 
value. In order to analyze the effect of the Hartree term 
we study also the functional Eq = E x + E z , that is the 
classical version of the quantum Heisenberg Hamiltonian 
studied by Henelius et al. jl3| In that work the Hartree 
term is not taken into account. 

A comment about the significance of is in order. 
In two dimensions the exchange term and the topologi- 
cal charge are scale invariants and do not depend on a^. 
The Zeeman energy increases quadratically with the lat- 
tice parameter and Vij is inversely proportional to a^. 
An increment in is similar to an increment of the Zee- 
man strength and the lattice parameter acts as a low 
energy cutoff that controls the dynamics of the classical 
Heisenberg model. We obtain the value of by fitting 
the magnetization obtained from the functional Eq to the 
magnetization obtained from quantum MC simulations. 
p3| In this way we obtain a temperature dependent a^. 
In the T — > 00 limit should be a\ . 

Now we describe briefly the MC procedure used for 
obtaining M(T). The MC simulations were performed 
by using the techniques due to Metropolis et al. C^] we 
study a cluster N x N with periodic boundary conditions 
(PBC). The use of PBC diminish the finite size effects. 
For studying MiT) at v = 1 we consider as a start- 
ing configuration a completely ordered ferromagnet, i.e. 
f2i = (0, 0, 1) for all sites i. In this way the initial Q in the 
system is zero. The sites to be considered for a change 
in the spin orientation are randomly chosen, avoiding ar- 
tificial correlations that could distort the results. Once 
a site is selected for a spin reorientation, we perform the 
following operations: i) a small change in the direction 
of fij, ii)calculation of the change induced in the topo- 
logical charge by the spin reorientation. Only changes 
that do not modify Q are accepted. iii)evaluation of the 
energy variation AE. iv) acceptance or not acceptance 
of the new spin direction. If AE < the change is ac- 
cepted. If AE > the change is accepted depending if 
e -AE/kT exceec j s a random number. We realize a large 
number of MC steps until reach the equilibrium situa- 
tion and perform the average of the different statistical 
properties: the magnetization M(T), the different contri- 
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butions to the total energy per electron, < E x > , < E z > 
and < E c > and a measure of the charge fluctuations, 

In Fig.l we plot the magnetization as a function of 
temperature as obtained from the classical MC simula- 
tion for the functional Eq = E x + E z with different value 
of the lattice parameter. The results correspond to a 
2DEG of zero layer thickness, p — 0.0249e 2 /i? and a Zee- 
man coupling t = 0.008e 2 /£. The calculations are per- 
formed in a cluster 20x20. We have checked that for 
this cluster size and PBC, the results are free of size ef- 
fects. In the same figure we plot M(T) as obtained from 
quantum MC simulations. |l3). This calculation does not 
include the Hartree term and it is the quantum version 
of the functional Eq. By comparing the classical and 
the quantum results we can estimate ar, in the different 
temperature ranges. We are interested in temperatures 
in the range 0.075e 2 £ < T < 0.2e 2 /£ At lower tem- 
peratures the quantum effects are very important and 
at higher temperatures the effective functional Eq.(l) it 
is not longer valid. In order to describe the dynamic 
in this range of temperatures values of aj, in the range 
2a± < aj, < 2.5ai are necessary. For simplicity, in the 
calculation we use a constant value of the lattice param- 
eter. We use aL—2.5ai, which is the value of ol in the 
range 2a\ < ol < 2.5ai where the effects of the Hartree 
interaction are weaker. 

In Fig. 2 we plot M(T) as obtained from the classical 
MC simulation for the full functional E — E x + E z + E c 
and for the functional without Hartree energy Eq = 
E x + E z . the results corresponds to a lattice parame- 
ter ar, — 2.5ai. At low temperatures, T < 0.05e 2 /£, 
the spins are not very disordered, the charge modula- 
tion is very weak and therefore the Hartree term has a 
small effect on M(T). For higher temperatures the spin 
fluctuations and consequently the charge fluctuations are 
stronger and the Hartree term becomes a important con- 
tribution to the internal energy. We find that at interme- 
diate temperatures, in the range 0.05e 2 /£ < T < 0.15e 2 £ 
the magnetization obtained with the full functional E 
is near thirty per cent higher than the obtained with- 
out the Hartree term. At even higher temperatures, 
T > 0.15e 2 /£, the spin disorder is very large and the 
magnetization calculated with or without Hartree term 
is small, although M(T) obtained with E is always higher 
than the obtained with Eq. 

As commented above the quantum MC calculation de- 
scribe correctly the spin density waves and at low tem- 
peratures it should give an appropriate M(T) for the 
v = 1 Hall ferromagnet. However we expect that for in- 
termediate temperatures the Hartree energy term would 
modify the quantum MC data in a similar amount that 
it modifies the classical results. In order to understand 



the experimental results at intermediate temperatures it 
is necessary to take into account the charge fluctuations 
induced by the temperature. 

In figure 3 we plot the different contributions to the to- 
tal energy per electron as a function of the temperature. 
The parameters are the same than the used in figure 2. 
At very high temperatures the spin are completely ran- 
dom and the exchange and Zccman energies tend to their 
fully disorder values 2p/o 2 and t respectively. Note that 
intermediate temperatures the Hartree energy is near one 
third the Zeeman energy. For smaller values of the Zee- 
man coupling and high temperatures we have found that 
the Hartree energy can be the more important contribu- 
tion to the internal energy. 

Figure 4 shows the variation of the charge fluctuations, 
Sq as a function of the temperature. Charge fluctuations 
cost Hartree energy and they are weaker when the full 
functional is considered. Observe that at intermediate 
temperatures 8q is of the order of 0.1, that is ~ 10% of 
the charge in each cell. 

In closing, we have studied the effect that the Hartree 
energy term has on the temperature dependence of a 
2DEG at v = 1. We find that at intermediate tempera- 
tures the spin fluctuations are weakened by the Hartree 
energy and the magnetization is near thirty per cent big- 
ger than the obtained by neglecting the Hartree energy 
term. At intermediate temperatures the Hartree energy 
is an important contribution to the internal energy of the 
Hall ferromagnet. 

We thank A.H.MacDonald, C.Tejedor, L. Martin- 
Moreno and J.J.Palacios for useful discussions. This 
work was supported by the CICyT of Spain under Con- 
tract No. PB96-0085 and by the Fundacion Ramon Are- 
ces. 



[1] A.H.MacDonald, H.A.Fertig and L.Brey, Phys.Rev.Lett. 

76, 2153 (1996). 
[2] For a review of quantum Hall ferromagnets see, 

S.M.Girvin and A.H.MacDonald in Novel Quantum 

Liquids in Low-Dimensional Semiconductor Structures, 

S.Sas Sarma and A.Pinczuk (eds.) (Wiley, NY 1996). 
[3] S.E.Barret et al. Phys.Rev.Lett. 72, 1368 (1994); 

R.Tycko et al. Science 268, 1460 (1995); S.E.Barret et 

al. Phys.Rev.Lett. 74, 5112 (1995). 
[4] S.L.Sondhi et al. Phys.Rev.B 47, 16419 (1993); 

H.A.Fertig et al. Phys.Rev.B 50, 11018 (1994); 

H.A.Fertig et al. Phys.Rev.B 55, 10671 (1997); L.Brey 

et al. Phys.Rev.Lett. 75, 2562 (1995). 
[5] E.H.Aifer et al. Phys.Rev.Lett. 76, 680 (1996); 

A.Schmeller et al. Phys.Rev.Lett. 75, 4290 (1996); 

V.Bayot et al. Phys.Rev.Lett. bf 79 1718 (1997). 

D.Leadly et al. Phys.Rev.Lett. 79, 4246 (1997). 
[6] M.J.Manfra et al. Phys.Rev.B 54 17327 (1996). 



3 



[7] 
[8] 

[9; 

[10 

in 

[12 
[13 
[14 

[is; 

[16 
[17 



N.Read and S.Sachdev, Phys.Rev.Lett. 75, 3509 (1995). 
C.Timm et al. Phys.Rev.B 58, 1464 (1998). 
M.Kasner et al. Phys.Rev.Lett. 76, 3204 (1996). 
R.Haussman Ph ys.Rev.B 56 9684 ( 1997). 



M.Kasner et al. cond-mat/9808186 



T.Chakraborty e t al. Europhys.Lett . 38, 141 (1997) 



P.Henelius et al. |cond-mat/9903148 



R.Rajaraman, Solitons and Instantons, (Nort Holland, 
Amsterdam 1982). 

M.Abolfath et al. Phys.Rev.B 56, 6795 (1997). 
B.Berg et al. Nucl.Phys. B 190 412 (1981). 
N.Metropolis et al. J.Chem.Phys. 21, 108 (1953). 
1.0- 




0.05 0.1 0.15 

KJ (e 2 A) 



0.2 



FIG. 1. M(T) as obtained by using the functional Eo with 
different lattice parameters. We also plot the quantum MC 
(QMC) results. [13] The results correspond to a 2DEG of zero 
layer thickness and a Zeeman coupling t = 0.008e 2 /£ 
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FIG. 3. Different contributions to the internal energies 
per electron. 
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FIG. 4. Variation of Sq, equation (8), as a function of the 
temperature for the full functional E and for the functional 
without Hartree energy Eo. 



FIG. 2. M(T) as obtained by using the functional E and 
as obtained by neglecting the Hartree term. The results cor- 
respond to a lattice parameter ah = 2.5\/2ty£. 
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